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Abstract: We investigate the generating functions of multi-colored discrete disks with 
non-homogenous boundary condition in the context of the Hermitian multi-matrix model where the 
matrices are coupled in an open chain. We show that the study of the spectral curve of the matrix 
model allows to solve a set loop equations to get a recursive formula computing mixed trace 
correlation functions to leading order in the large matrix limit. 



1 Introduction 

The theory of random matrices has many ramifications in different fields of mathematics and physics. 
In the recent years, following the work of Eynard [1 Ij , matrix models have seen significant advances 
especially through their applications to enumerative geometry. The enumeration of surfaces can indeed 
be addressed by the evaluation of the result of a saddle point approximation on integrals defined over 
a set of random matrices as the size of the matrices gets very large. To be more precise, the result 
of such an approximation can be seen as the generating function of a set of fat graphs composed of 
ribbons glued together along vertices [J. This set of fat graphs is in bijection with the set of discrete 
surfaces, or maps as denoted by combinatorists, i.e. surfaces composed of polygons glued by their 
edges. The matrix integrals can thus be seen, in this large matrix limit, as the generating function of 
a set of maps with a weight depending on the probability measure put on the set of matrices under 
study [UirniH^. Such objects are of prime interest not only to mathematicians but also to physicists 
for investigating the behavior of statistical systems on random surfaces (or lattices) as well as for 
understanding quantization of two dimensional gravity or topological string theories from a discrete 
point of view. Discrete surfaces also provide a very good toy model for developing some intuition to 
be applied later to string and gauge theories. 

From this perspective, a lot of progress have been made by Eynard et al. leading to powerful 
universal technics for the computation of generating functions of maps of different types. Among the 
results obtained, one finds: 
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• The enumeration of maps of arbitrary topology, i.e. arbitrary genus and number of boundaries, 
has been achieved in [U [Sj [11] ; 



• The enumeration of bicolored maps of arbitrary genus and arbitrary boundary conditions has 
been completed in the series of papers [HI [HI [HI [TSl [HI [TS] ; 

• The enumeration of bicolored maps drawn on non-oricntable surfaces has also been solved in 
[51 [ini [3 [Hl|3] with homogenous boundary conditions; 

• The enumeration of maps whose faces are colored with more than two colors has been partly 
solved in |131 123) . In particular, the generating functions of maps with homogenous boundary 
conditions are known whatever their genus is. 

One big outcome of these works is the universality of the solution found for all these problems. 
Indeed, for each of the aforementioned enumerative problems, the generating function of maps with 
homogenous boundary conditions are given by a unique universal inductive formula on the Euler char- 
acteristic of the surfaces enumerated. Even better is the fact that this inductive formula, sometimes 
called topological recursion, goes beyond the field of random matrices and enumeration of maps. The 
topological recursion |201 122j seems to be a universal solution to different problems of enumerative 
geometry consisting in the computation of the volume of the moduli space of Riemann surfaces with 
respect to different measures opening the way for new insights in mathematics and topological string 
theories even when this measure does not localize on discrete surfaces. 

The topological recursion was inspired by the computation of generating functions of open surfaces 
with homogenous boundary conditions in the different matrix models encountered. The next step in 
fully understanding the enumeration of maps in general is to understand the action of boundary op- 
erator, that is compute the generating function of surfaces with boundaries along which the boundary 
condition changes. This has been understood in the hermitian two matrix model case where it has 
been possible to compute the generating function of surfaces composed of polygons of two different 
colors with arbitrary boundary conditions [19] . The result points towards a generalization of the topo- 
logical recursion for these new observables. Going further requires the study of the so-called chain of 
hermitian matrix model which generates maps colored with an arbitrary number of colors. Even if 
the homogenous boundary condition case is understood [55], there are only a few result concerning 
changes of boundary conditions [13] . 

In these notes, we present a recursive method, similar to the topological recursion, which allows 
to comp ute t he generating function of discs with mixed boundary conditions. In theorem |4.1| and 



can be computed to leading order in the large matrix limit in terms of the spectral curve of the model. 
The formula takes a very nice form which seems to follow from the all possible degeneracies of the 
surfaces generated. From the two matrix model experience, it is probably the building block for all the 
other generating functions and probably one of the most important objects of the theory. We expect 
to be able to use the techniques developed in this paper to compute any observable of this model. 

The computation of these disk amplitudes has thus many interesting features. First of all, it gives 
an efhcient method for computing the generating functions of discs composed of colored polygons 
with mixed-boundary conditions. In addition, we are convinced that the result as well as the methods 
used in its derivation are generalizable to the computation of the generating function of surfaces 
of arbitrary topology. From this perspective, these notes should provide the building block for the 
computation of all observable of the chain of matrices model. This result is also important for better 
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understanding conformal field theories. Indeed, when going to a particular limit in the parameter 
space of matrix models, one reaches critical regimes described by minimal models, a particular type 
of rational conformal field theories. The knowledge of boundary operators in the matrix model setup 
thus also gives access to boundary operator in such rational conformal field theories. Finally, the form 
of the solution points towards a generalization outside of the matrix model's setup generalizing the 
topological recursion formalism leading to the knowledge of new open amplitudes in topological string 
theories or may be new relative Gromov-Witten invariants. 

This paper is organized as follows: 

• In section 2, we introduce the model as well as the notations required in the following and 
present the result of these notes; 

• In section 3, we remind how algebraic geometry appears in this context through the so-called 
loop equations and introduce all the notations and basic knowledge further needed. 

• In section 4, we derive and solve a closed set of loop equation satisfied by the leading order 
mixed correlations functions. 

• In section 5, we specialize this result to the chain of 3 matrices and give explicit examples. 

• Section 6 is devoted to the conclusion and a short discussion on the perspectives. 



2 The open chain of matrices 
2.1 The model 

Wc want to generate discrete surfaces, i.e. surfaces composed of polygons glued by their edges. In 
addition, we want to include some additional information on these random surfaces such as a spin 
structure or some matter. For this purpose, each polygon carries one color labeled by an index ranging 
from 1 to M. 

Generating functions for such discrete colored surfaces can be obtained as correlation functions 
of a multi-matrix model. The latter studies the distribution of matrix elements with respect to a 
probability measure rf/i on the set (i?Ar)'^ of AA-uple of hermitian matrices of respective size N x N 
defined as 

(2-1) 



/ N 

diJ.{Mi,...,Mj^) :=dMidM2...dMj^exp Tr 



T 



■ AT M 
^ VkiMk) - ^ Ck-i,kMk-iMk 

.fe=l k=2 



where Mi is a hermitian matrix of size N, dMi is the product of Lebesgue measures of the real 
components of Mj and 

dk + l (fe) 



i=2 



are polynomial potentials of respective degrees (ifc + 1. The partition function of this model is defined 



as 



Z:= / dii{M^,...,M^). (2-3) 

In this paper, we interpret this random matrix integral as a generating function of discrete surfaces. 
This interpretation is exact if one considers a formal version of this integral, which is defined as 
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the result of a saddle point approximation of this on^ In the following, all the matrix integrals 
encountered are formal in this sense. 

In a nutshell, the combinatorial interpretation of such formal matrix integral is obtained by ex- 
panding the non quadratic part of the action 



fc=l k=2 



(2-4) 



around the extremum {Mi — 0}^j^ and represent the result by Feynman graphs whose edges are 
ribbons which cannot be twisted. In the expansion of the action around this saddle, each factor 
Tr Mj' is represented by a fc-valent vertex of color i whose legs are ribbons. This peculiarity of the 
edges puts an ordering on the legs around every vertex. One can thus replace every vertex by a face 
(or polygon) whose number of edges equals the valence of the initial vertex. Then, each graph is dual 
to a map, or discrete surface. 

Hence, the partition function of the matrix model is defined as the generating function of closed 
surfaces made of polygons colored with J\f different colors and glued together following the prescrip- 
tions: 



• the polygons of color k have at most dk + 1 edges ; 

• two polygons of respective colors i and j can be glued by their edges only. 



In this representation, the partition function reads 



oo 



oo 
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where Mg(w) is the set of orientable genus g maps with v vertices, Aut(S') is the number of automor- 
phism of S, rii^kiS) the number of faces of color k with i edges in S, nij{S) the number of edges 
between a face of color i and a face of color j in S and [C'""'^] the inverse of the matrix of the bilinear 
terms of the action 



C :-- 



The free energy 



V 
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(2-6) 



T := InZ 



is the generating function of the connected surfaces contributing to Z. 

For any function / : (Hpf)'^ — > C, one defines a corresponding correlation function by 



(/(Mi,...,AfAA)) := 



1 



fiMi,...,M:^)df,{M, 



(2-7) 



(2-8) 



The correlation functions of this model can also be interpreted as the generating functions of open 
surfaces built with the same prescription and the same measure as the closed ones, i.e. surfaces with 



^This formal integral is defined as a formal scries in T and docs not need to converge. More precisions about this 
definition can be found in 1141 
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Figure 1: Example of map contributing to ( Tr MJ* Ml M|) where the type 1 condition (resp. type 
2 and type 3) is depicted in yellow (resp. red and blue). It is important recalling that the boundary 
condition corresponds to the color of the outer face which has been removed and not the color of the 
inner polygons touching the boundary. The dual fat graph is depicted in dotted lines. 



boundaries and prescribed boundary conditions on it. Let us briefly recall how to define generating 
functions of such surfaces. Consider a correlation function of the form: 

( Tr \{M,, (M,, . . . {A'U, 1 \ := 1 /■ Tr [(M,, ^ (M.. f ■ ■ ■ (M^, ]d^i{M,,..., M^) 
\ L J/ 2 J^H^)^ L J 

(2-9) 

where (ii, «2, . . . , «fc) G [liA/"]*^ and (ji,...,jfc) £ (N*)' are arbitrary positive integers. The factor 
Tr (Mi-^Y^ (^22)"'^ ■ • ■ (-^ifc)""' in the integrand constrains all the generated graphs to contain one 

k 

vertex with ii colored legs organized so that the sequence of colors coincide with the sequence 

1=1 

of matrices in (M^J"'^ {Mi^Y^ ■ ■ ■ {Mi^.Y'' ■ Thus, this correlation function is the generating function 
of all possible fat graphs that can be glued to this vertex, i.e. all fat graphs with one boundary 
constrained by the condition that it must be glued to this vertex. In the dual representation where 
k-valent vertices are replaced by k-gones of the same color, this function generates all surfaces with 
one boundary constrained by the same type of gluing condition. 

In order to deal with all the possible exponents ji at once one introduces the resolvents 



00 

^31 + 1 32 + 1 + 1 



n,^.=i (2-10) 
■= i/xr 1 1 i \ 

as formal power series in its parameters xi, 0:2, . . . , x^. It is very important to note that these resolvents 
are asymptotic series defined only when there arguments are large!. Indeed, we will show that these 
correlation functions are multivalued functions of their variables x in the complex plan^ but can be 
promoted to monovalued functions on an Riemann surface. However, one must always keep in mind 
that the physical quantities are the coefficients of the expansion of these correlation functions around 
a very particular point. 

Using the observation of 't Hooft for QCD one can also select the genus of the generated 
surfaces. In order to do so, one can show that these correlation functions are series in and that 
the order of is genus of the generated surface, i.e. one can write: 

00 ^ 

W,,,,,, (Zi, Z2,...,Zk)=J2 -]^^^u^2,...,^, (^1' ^2, • • ■ , ^fc) (2-11) 

9=0 



^This property can be seen as the result of the finiteness of the radius of convergency of this series. 
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where W^^l^ (zi, Z2, • • ■ , ^fc) is the generating function of genus g surfaces with one boundary with 
the boundary condition induced by Wi^^i^,,,,^if.(zi, Z2, ■ ■ ■ , Zk)- In this paper, we will be interested in 
the large N limit of these correlation functions, that is to say, the generating functions of discs with 
boundary operators of the form W^^\^ (zi, Z2, . • . , Zfc). 

2.2 Main result: disk amplitudes 

The main result of this paper is the computation of the mixed disk amplitudes with the use of a 
recursion relation explained in th|4.1| 



Wi,j,j+i,...,7v(pi,Pj, • • ■ ,P7V) = Res K:■j{pl,q■j,pJ)Wl^■j+l,...,^riq,pJ+l,. . . ,p^r) (2-12) 

where the ingredients are explained in the following of the paper. In particular, this formula makes a 
heavy use of the different sheeted structures of the spectral curve of the model under study. 

Another important feature of these notes are the explicit computation of these amplitudes as well 
as their expansion in the 3-matrix model. In section [5] we apply the spectral curve approach as well 
as th 4.1 to explicitly compute the number of tri-colored discs with given mixed boundary conditions. 

2.3 More notations 

In this section, we present the notations used in the following of these notes. We used the notations 
of |23j as long as it was appropriate so that the reader interested in both aspects of this matrix model 
can go from one paper to the other without being confused. 

First of all, let us remind the correlation functions to be computed: 

T / 1 1 1 \ 

iV \ xi~ Mi^ X2 - M,.^ Xk - Mi^ I 

We also define the generating function of surfaces with two boundaries by 
W'i-.7(a;i,a;2,...,Xfe;2/i,...,?;;) := ( Tr ( - — ^^ri- • ■ • - — ) 



(2-14) 



^Xi-Afji Xk-Mi^J Vyi-Mjj yi-Mj^ 
where the index c stands for the connected part and 

i = (ii,Z2, . . . ,ifc) , j = {ji,j2,---,ji) (2-15) 

are sequences of integers in [IjA/"]. 

Following |13j . let us now define a set of polynomials fij{xi, . . . ,Xj) by induction: 

( := if J<i-1 

fii—l '.— 1 

y (^\ — Yli^ ■ (2-16) 

For 1 < fc < Z < Af, let us then define 
Pk,i;ji,...,j^ixi, ■ . .,xj^;yi, . . ■,ym) ■= Pol fk,iixk, ■ ■ ■ , Xl)Wl^,„^^r■J^J^^„„j^{xl, . . .,xj^;yi, . . .,ym) 

Xk,---:Xl 

(2-17) 

which is a polynomial in the variables x^, i = k, . . . ,1. 
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Given a pair of complex numbers {xi, X2), we define a set of polynomial functions {xi{xi, X2)}'^^ 
by the recursion relation 

Ci^,+iXi+i := Vl{xi) - Cj_ii;j_i (2-18) 
for alH = 2, . . . ,J\f. The case i = 2 is slightly different and we define 

Y{x,) ■.^V;ix,)-W,{xi). (2-19) 

It should be noted that the latter function has a expansion. 

Finally, it is convenient to denote any observable whose variables satisfy the constraints Eq. (2-181 

by: 

0{xi,X2) ■■= 0{xi,X2,Xz{xi,X2),Xi{xi,X2), ■ ■ . , x^f {xi, X2)) ■ (2-20) 

All the observables of the model admit a topological expansion [2^, i.e. they can be written as a 
formal power series in 

oo 

O ^^N^^SQia)^ (2_21) 

3=0 

3 Master loop equation and spectral curve 

Before going into the explicit computation of correlation functions, let us recall the algebraic structure 
underlying the hermitian multi- matrix model studied by Eynard in 13 . We will particularly focus 
on the sheeted structure of the spectral curve since it is the most important ingredient needed in the 
following. The interested reader can find a longer and more complete presentation of this structure 
in |13J or more recently in |23j . 

3.1 Loop equations 

The partition function and the correlation functions, seen as generating functions of random surfaces, 
are constrained by an infinite set of equations prescribing how the weight of a surface is changed when 
one remove one edge from it. These equations where first obtained by Tutte [571 HH] for counting 
triangulated surfaces and later introduced under the name of loop equation in the context of matrix 
models by Migdal [24]. These equations proved to be some very powerful tools since they allowed the 
computation of many observables in different types of hermitian matrix models through the arising of 
algebraic geometry. 

Practically, these equations can be obtained by writing the invariance of the integral in Eq. ( 2-3 ) 
under particular infinitesimal change of variables: 



■ Tr 



TV AT 

Vk{Mk) - Ck-i,kMk-iMk 

k=l k=2 



Z = J dMi . . . dMj . . . dMj^e 

= / dMi . . . dMj . . . dM^r J{Mj,Mj) X 

-fTr Y'^k{Mk) + Vj{M,)- Y Cfe-i,fcMfc_iMfc - Cj_ijAfj_iMj - Cjj+iA?jMj+i 

(3-1) 

where Mj := Mj + eSMj is an infinitesimal deformation of Mj. The loop equation just corresponds 
to saying that the variation of the action is compensated by the Jacobian of this change of variable 
to first order in e. 
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The set of loop equations corresponding to the changes of variables 

Mfc Mk + e ^- ^- . . . ^- ^- (3-2) 

Xk+i - A/fc+i Xk+2 - Mk+2 - M^f xi - Ml 

supplemented with some algebra allows to show that 
Theorem 3.1 For < k < M, the correlation functions satisfy 

{xk - Xk{xi)) T4^i,fe,...,7v(a;i, Xfc, . . . , Xf^) = W^i,fc+i,...,7v(a;i, x^+i, . . . , x^f) 

-Pl.k-lixi,X2, ■ ■ ■ ,Xk-l,Xk, . . . , Xf^) 

T 

+ J^P2,k-l-l[Xl,X2, . . . ,Xk~l,Xk, ■ ■ ■ ,Xjs/) 

(3-3) 



where 



Pr,k{xi, . . . ,xaa) := ^ Pol^^ fi,k{xi,. ■ ■ ,Xk-i)Wi^,,,x{xi, . . . ,X7v) (3-4) 



is a polynomial in the variables Xj, j — i,. . . ,k and the polynomials fi^k{xi, . . . ,Xk) are defined by 
induction through 

J. I N Vki^k)fi,k~l — Ck-l.kXk~lXkfi.k~2 /„ 

fi^k\Xi,...,Xk-i) = ■ ■ . (3-5) 

Cfe,fe+1 

For k = N , one has 

(Kv-(^aa) - ^TV-i) Wi{xi) = Pi,m{xi,x2, ■ . ■,x^) + j^P2Mai^i'^2, ■ ■ ■,xm;xi). (3-6) 

This theorem was proved in 1131 123j . For completeness, it is derived for the present model with 
the notations in hand in appendix A. 

3.2 Leading order and master loop equation 



By plugging Y{x) ~ Vl{x) — TWi(x) into Eq. (3-3l, the leading order in the expansion of this 
equation reads 

E{xi,y{xi))=0. (3-7) 

where 

2/(x)=y(")(x) (3-8) 
is the leading order of the expansion of Y{x) and 

E{xi,X2, . ■ .,x^r) = {V({xi) - x2)iV^r{x^f) - xw-i) - Pi%{xi, ■ ■ -.x^f) (3-9) 

is a polynomial in all its variables. 

Consider Eq. p-7[ ) as an equation with unknown xi and y{xi). On can see that this is a polynomial 

AT J\f 

equation of degree JJ^i in xi and 1 + Wfii m y{xi). Thus, it defines an algebraic curve £{x,y) = 

i=l i=2 

in C\ 
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3.3 Description of the spectral curve 

Let us now consider the algebraic equation 

Eix,y) = 0. (3-10) 

It defines an algebraic curve in which encodes most of the combinatorial properties of the matrix 
model: it is the spectral curve associated to the model. 

In this section, we study the properties the algebraic curve E{x, t/) = as a compact Riemann 
surface C equipped with two meromorphic functions x and y such that 

ypeC, E{xip),y{p))=0. (3-11) 

A point p G C is equivalent to a pair {x(p), y{p)) G satisfying E{x{p), y{p)) = 0. 

Since, given a pair (x,y) G C^, all the Xi{x,y) are unambiguously fixed, one can see them as 
monovalued meromorphic functions Zi on C such that: 

Vpe£, Zi{p) = Xi{x{p),y{p)). (3-12) 

By definition, one has 

x{p) = zi{p) , y{p) = Z2{p). (3-13) 
Using these functions, one can write: 

Vpe£, Eiziip),Z2ip),Z3ip),...,z^{p))^0 (3-14) 

and one can choose to fix any of the functions Zi{p) as parameter to study the curve. Depending on 
which parameter we chose we obtain different descriptions that we can now characterize. 

3.3.1 Sheeted structure and points at infinity 

Let us first study this curve in terms of the "physical" variable x := Zi. Since E{x, y) is a polynomial 
in y of degree 1 + y{x) is a 1 + J^di-valued functions of x. This means that, for a generic value 

i=2 i=2 
N 

of X, there exist 1 -|- JJc^i points p> on C such that: 

i=2 

N 

= 0,...,JJdi , zi{p') = Zi{p). (3-15) 

i=2 

Hence C can be seen as 1 -|- si copies of the Riemann sphere, called 2;i-sheets, glued together by 
cuts in such a way that x{p) = zi{p) is injective in each sheet. This sheeted structure reflects the 
multi-valuedness of the resolvent W^"^ {x) as a function of x. Indeed, once the value of this correlation 
function is fixed, y{p) = Z2{p) is fixed and all the Zi{p) = Xi{x{p),y{p)ys are determined. Thus, this 
fixes one unique point in £. 

How can we extract the right value of y{x) giving access to the generating functions of random 
surfaces? In order to distinguish the different a;-sheets, one can look at the pre-images of oo by x{p), 
i.e. the points oo-' e jC such that a;(oo-') = oo. 
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From its definition as an asymptotic series, one knows that the physical sohition of the algebraic 
equation should satisfjQ 

y{p) ^.(p)^oo V'{x{p)) + ^ + 0(x-2(p)). (3-16) 
' x(p) 

It means that there exists at least one pre- image of infinity, denoted oo^ e £, such that 



1 

zi{p) 



Z2{p) ^p^ooi Vi{z^{p)) + + 0{z^\p)) (3-17) 



and more generally 

Vfc = 2,...,AA, Zfe(p)^p_ooi O(zJ^) (3-18) 

where 

One could symmetrically choose zjv as a variable and use the physical condition that 

ZM-l{p) V;^{zm{p)) + —r-, + 0{z^\p)). (3-20) 

ZMKP) 

This implies that there exist at least one pre- image cxd-'^ of infinity by x such that: 

zi{p) -p-.oo^ 0{z^{pY') (3-21) 

where 

rk = X{i^k + l^d,. (3-22) 

This means that ri different zi-sheets merge at oo-'^. Since ri + 1 is exactly the number of zi-sheets, 
there exist no other pre-image of infinity by x = zi and the fiber above infinity reads: a;~^(oo) — 
{ooi,oo^}. 

One can summarize the structure of E in terms of zi as followf[^ 

• £ is composed of 1 -I- ri zi-sheets; 

• zi[p) has two poles in C: one simple pole noted oo^ and one pole of degree ri noted oo-'^; 

• oo^ belongs to only one sheet called the physical sheet since it corresponds to physical solutions. 
All other sheets merge at oo^ . 

The unique zi-sheet containing cxi^ is only referred to as the physical sheet. Indeed, in order to 
recover the physically meaningful quantities Tr , one has to expand y{z) around z — > oo^: 

viz) -.^oo^ Viixiz)) + Y: ^ { Tr M^) . (3-23) 

k=0 ^ ' 

The expansion of y{z) around oo^ doesn't have such a simple combinatorial interpretation. 

One can also fix Zk{p) and look at the Zfe-sheeted structure of C. With the same arguments, one 
can see that: 

• There exists + Sk Zfe-sheets where Zk{p) is injective; 



■^We call this solution physical because it is the branch whose expansion in x gives rise to the generating function of 
surfaces. 

^An example is depicted in fig[2] 
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Figure 2: Example of sheeted structure of L where one has represented a genus spectral curve as 
the Riemann sphere, the red dot represents the pole oo^ while the black dot is the other pole oo-'^. 
The different sheets are represented by the connected components separated by continuous lines. In 
this case, ri = 4 of them join at oo-'^ and only one contains oo^. 



• z-k(j)) has two poles. One pole of degree Sfe in oo"'^ where Sfc z^-sheets meet and one pole of degree 
Tfe in oo''^ where the other Zfe-sheets merge. 

• One notes p^^'^ (resp. p^^''^) for j — l,...,Sfc (resp. j = l,...,r/j) the pre-images of Zk{p) 
belonging to the different sheets which merge at oo^ (resp. oo-'^). 

In particular, one can see that there exists a neighborhood of oo^ where zi is injective in the 
Zfc-sheets. This means that: 

Zkip^'-'') = z,{q+^-'^) ^ x{p+^-'^) = x{q+^-'^) (3-24) 
in some neighborhood of oo^. 



4 Correlation functions 
4.1 State of the art 



Let us first recall the results of [131 123] . 

In [23', a very efficient inductive method has been found for computing the whole topological 
expansion of the free energy as well as the correlation functions involving only the first matrix of the 

chain. This means that the free energy as well as the correlation functions of the type ( JJ^ Tr ) 

are known to any order in the large N expansion. 

In [13], some other observable were computed. First of all, the leading order in the large N 
expansion of any two point correlation function of the type ^ Tr Tr ^ was computed 

and expressed in terms of a fundamental one form on C x £ often called the Bergman kernel. More 
important to us is the derivation of the simplest disk amplitude with mixed boundary conditions, that 
is the large N limit of the correlation function ( Tr ^ ^ 



4.2 Mixed correlation functions 

With this description of the algebraic curve in hand, we are ready to compute the correlation functions. 
Let us first promote them to mono-valued differentials Wij^,...i^ (pi, . . . ,pk) on C*' defined by: 

y^ii,...tk (Pi, ■■■,Pk) ^ ■ • ■ , {Pk)) dz^^ipi), dz,^{pk). (4-1) 

We use the same notation for any other function considered in the previous section; the curly letters 
represent differentials on the spectral curve. 
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From now on, we will consider only the leading order of the correlation functions in the 
expansion. We thus abusively denote 

Wii,j2,...,,^^(xi,a;2, . . . , Xk) := W^^l^^_^^^{xi, X2, .■■,Xk) (4-2) 

and 

Pk,i{xi, . . ■,xj^) := Pfc°^(a;i, . . .,xj^) (4-3) 
the leading orders of the observables we are studying. Let us also remind that 

yix,) := r(°)(xi) = Vlixi) - Wi°\xi). (4-4) 

In terms of differentials on the spectral curve, to leading order in the expansion, the loop 

equation Eq. (3-3) can then be written: 

i^jiPj) ~ Zj{pl))WlJ.J+l,...,^f{pl,pJ, . . . ,pj~j-) = Wi,j+i,...,A/-(pi,_Pj+i, . . . ,p^)dz-i{pj) 

-'Pi.j-i{pi,z2{pi), Zj_i(pi),pj, . . . ,p^f) 

(4-5) 

where 

'Pl,3~l{Pl,Z2{pi),---,Zj-l{pi),Pj,---,PAr) „ / / N f \ f \ r W a\ 

— = Pij-i(2i(pi), . . . , Zj^i(pi),zj{pj), z^{pf^)) (4-6) 

is a polynomial in zi(pi) of degree Vj — 1. 

The main result of these notes is that one can solve this equation by induction on j. Indeed, 
Eq. (4-5) expresses >Vijj+i,...,AA(pi,Pi, • ■ • tPm) in terms of 'Wi^j+i^...^j^(j)i,Pj, . . . ,p^f) provided that 
one knows the polynomial Vij-i{pi, Z2(pi), . . . , Zj__i(pi),pj, . . . ,pj\f). The polynomiality of the later 
actually allows one to compute it by using only Eq. (4-5) and get 

Theorem 4.1 Recursion relation for the disk amplitudes. 

The mixed disk amplitudes satisfy the recursion relation 



yVi,j,j+i,....Af{pi,Pj,---,PAf) = R-es ''"^ 



[zi{q) - zi{pi)] [zj{pj) - Zj{pi)] Y[ [zi{q) - zi{p^'''^ 

(4-7) 

Since the initial value 



m,A.(Pi,P^) = E / .X/ ^ (4-8) 

fe=l \Zjk — ^h] I + 1 — ^k+1 



Zk ^ Zk (pt'^izi)j Zk = Zk {p^^ {z^/)^ (4-9) 



was computed in [13 , this theorem allows to compute Wijj+i,....A^(pi,Pj, . . . ,paa) for any j 
2, . . . , A/" — 1 and in particular the complete mixed correlation function: 
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Corollary 4.1 The disk amplitudes read, for j = 2, . . . ,Af — 1: 

Wij,j+i,...,AAbi,Pj, ■ • ■ ,Paa) = Res Res ... Res TT /C^bi, <?a,Pa)Wi,AA(gAA-i,PAA) (4-10) 



1i + l 9AA-1 

a=] 



where one defines the recursion kernel 



9, r) = (4-11) 

[zi((Z) - ziip)] [z„(r) - zdp)] n - 21(7-+'=^")] 



k=l 

and the residue sign stands for 



= j, . . . ,7V - 1 , Res = Res -h V Res . (4-12) 

O-i (7.-— VTIi ' ^ _ + 



In turn, this gives access to any disk amplitude of the form 'Wi^^...s^ with 1 < Ji < «2 < ■ • ■ < «fc < 
W by considering the leading order of the expansion the complete disk amplitude as some pi 00^. 



4.3 Proof of the recursion relation 

For proving th |4.1[ one first needs a simple technical lemma. 

Lemma 4.1 {zj{pj) — Zj{pi))Wijj+i^,,,^j^{pi,Pj, . . . ,pj^) vanishes when pi — p'^^'\ for i — l,...,Sfc. 
proof: 

The proof follows from the combinatorial interpretation of yVijj+i_...^7v(pi,pj, . . . ,paa) when pi 
lies in the physical sheet. 

Indeed, for pi —J- 00^ (and not in any other patch of the spectral curve), the definition of 
Wijj+i,...,AA(pi,Pi, ■ ■ ■ ,P7V) reads 

^ i^jiPj) - ZjiPi)) Wl,ij+l,...,7v(Pl,Pi, ■ ■ ■ ,pm) = 

^ ^'^^'^ ' '^^''^^ { zM - z,ip,) - M, ■ • ■ .^(J-M^) ■ (4-13) 



One can thus write 



N{zj{pj) - Zj(pi))>Vij.j+i,...,Ar(pi,_Pj, . . . ,pm) 

1 1 1 

Tr 



(4-14) 



zi(pi) - Ml Zj+i{pj+i) - M,+i ■ ■ ■ z^{p^f) - M^f 

^ Zj{pi)~Mj 1 

zi(pi) - Ml Zj{pj) - Mj ' ■ ■ z_^f{p_^f) - Mtv 



Tr 



which vanishes when Zj(pj) = Zj{pi), i.e. pi belongs to the same -fiber as pj. Since this 
formula is valid only for pi in the Zi physical sheet, this additional constraints implies that 
(^j(Pj) - 2j(Pi))Wi,jj+i,...,7v(pi,Pi, ■ • • ,P7V) vanishes for pi = □ 



Proof of theorem 14.11 

Thus, one knows that the left hand side of Eq. (|4-5|), vanishes for pi = p^ 
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On the other hand, the second term of the right hand side, Pi is a polynomial in zi(pi) of 
degree Sj — 1. Thanks to the vanishing of the left hand side, one knows its value in the c?i . . . dj-i 
points, zi{p'^''^) for j = 1, . . . , Sj : 



(4-15) 

By Lagrange interpolation, one gets an explicit expression for this polynomial which one can plug 
back into Eq. ( 3-3 ) to get 

{zj{pj) ~ Zj{pi))Wij....x{pi,Pj, . . . ,pu) = yVij+i,...,A^-(pi,pj+i, • ■ • ,PN-) 



E 



fc=i 



Zl 



{pi)-^i{pT-')]l[['i{pp'^-^iip- 



(4-16) 



This last equation can be written under the form of a residue formula 



Wij,...x(Pi'-Pi' ■ • ■ 'PW) 



r 

Wlj + l,...,A/-('?,J3j + l, • • • ,PAA)n [^l(Pl) - Zi{p'^^''' 



Res 



k=l 



[21(9) - 21 (pi)] {zj{Pi) - Zj{j)x)\ n \zx(si) - zx{p]^'') 



(4-17) 



□ 



5 Example of application: the 3-matrix model 

In this section, we perform explicit computations in the case A/" = 3 resulting in the computation of 
three colored discs with mixed boundary conditions. 

5.1 Genus spectral curve 

For enumerating maps, the spectral curve has to be a rational curve, i.e. its genus has to be vanishing. 
Hence, there exists a global coordinate p such that the functions Zi are rational functions [2] 

Ti 

Vi = 1, . . . ,3, Zi(p) = ^ (Xi^kp^ (5-1) 



k— — s 



satisfying 



y/(zi(p)) - Z2(p) ^p^oo + O (p-2) (5_2) 

POLX,\ 

y^{.^M) - 22(P) ^p^o — + O (5-3) 

^3,-1 

F2'(;^2(p)) = zi(p) + 23(p), (5-4) 

where a\^\ = Q!3,-i ~ 7. The limits of the sums describing the parametrizations are obtained from 
the degrees of the potentials, and are given by s\ = r^i — 1, r\ = ^2^3, S2 = d\^ri = and S3 = did2. 
The solution branch giving rise to the spectral curve is the one where the coefficients 7 and ai^k are 
algebraic functions of T such that 7 — ^ as T — ?> 0. With this global parameterization ooi is located 
in p = 00 whereas 003 is p = 0. 
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Once the spectral curve, i.e. the collection of rational functions Zi(p), is known, one can study 
its different sheeted structures. The Zj-fiber over a point Zi is given by the set of solution p'^^''^{Zi) 
satisfying 

Zi {p^^'\Zi)) = Zi. (5-5) 
Among these solutions, oo and p~^''^{Z) — )• as Z oo. 

Remark 5.1 In the setup of enumerativc geometry, the exponent of T is the generating functions is the 
number of vertices in the maps considered. Thus, as long as one is interested in sufficiently small maps, the 
preceding equations need being solved only up to a given order in the T expansion. 
One can thus compute the fibers above a point Z through its small T expansion. 

5.2 The complete one loop function 
5.2.1 Loop equations and analytic solution 

Since the 3-matrix model is the simplest model which carries the whole complexity of the chain of 
matrices, let us use it as a toy model for deriving more carefully the loop equations involved in the 
computation of the complete mixed correlation function. 
Let us consider the following change of variables 

Ml Ml + e5Mi (5-6) 

where 

,5Mi = ?— . (5-7) 

In order to keep the path-integral invariant under such a change, we demand that the contributions 
from the Jacobian of the change of variables and the contributions from the variation of the action 
cancel each other. To order e, we have 

(^) ) = {{Tr{Vl{Mi)6Mi)} - {Tr{5M,M2))) (5-8) 
The left hand side can be written in terms of the connected correlators 

(^■-f*"''^" (57^)) = <^"(™'»= (^'' (jT^) + (^-<™-'^'' (^) 

(5-9) 

therefore we can write the left hand side in terms of correlators dependent on the points p, g, r e E as 

iV'Wi,2,3(p,«,r-)Wi(p) + Wi,2,3;i(p,g,r;p) (5-10) 

where 

zi{p) = zi , 22(g) = ^2 , zz{r) = zs. (5-11) 
On the other hand the right hand side can be expressed as 

- N {Tr{{Vi{Mi) - Vl{zi))SM,)) + N\V'{zi{p)) + 22(?))Wi,2,3(p,«,r-) - N^Wi,s{p,r). (5-12) 
By collecting the terms and considering only the order N° terms, we have 

{Z2{q) - Z2{p))Wi,2,3{P^q,r) = Wi,3(p,r) - Pol Fi'(zi(f>))Wi,2,3(p, g, r), (5-13) 

zi(p) 



where 



mMP,r) = Mfi^ _ Mfl^ (5.14) 

Z2 — Z2 \ z\ — Zi z-j, — Zj-i / 
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[C-'],3[C']33g3=2g3 



Figure 3: Maps contributing to Ti q^^- The type 1 (resp. 2 and 3) boundary condition is represented 
by the color white (resp. red and blue). The outer face is the marked face of type 1. 



is known from [13j, and where we used the notation: Z2 — V({zi) — Wi{zi), Z-i — V^iz^^ — 14^3(23), 
Z3 = V!,{Z2) - zi and = V^^Z^) - z^. 



This equation can be solved and its solution is given by th 4.1 under the form 

/ 

Wl,3(Q,P3 

fc=l 



Wi,2.3(Pl,P2,P3) 



Wl,3(Q,P3)n(^l(?'l)-^l(P^''')) 



:^i(o) - ^i(Pi))(^3(P3) - ^3(Pi)))n(^i(Q) - 



fe=l 



(5-15) 



5.2.2 Triangulations 



In this section we present some explicit computations for the three matrix model. We assume that all 
of the matrices are subject to cubic potentials of the form 



(5-16) 



Furthermore, for simplicity we assume g'^'^ — ~ 1 and g''2^ — 3 in order to obtain a convenient 
form for the propagator matrix: 



2 -1 1 
-1 1 -1 
1 -1 2 



(5-17) 



The conditions 5-2 can be used to express the ZiS as a -v/T-expansion whose coefficients are 
functions on the spectral curve that also depend on the coefficients of the potentials. 



n=l 

For example, the first few orders of zi{p) are given by 

z['^^gi^^;p)=-(p+l), 
which can be recognized as the gaussian contribution, 



(5-18) 



(5-19) 



(5-20) 



, 25^5^ , -325r)-6,«5r-2ffW.9f + 16(5^)2 
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-o 







[C'],,[c']jjc']^,gr4gi [C'],jc']^jc-']^,g2= [C'yc-'],3[c-]^^g3=2g3 



A 





[c-']„[c-']jc-'],^g,^2g, [c']jc-'yc']^^g,^ig^ [c']jc'yc-]^^g,.ig, 



Figure 4: Maps contributing to T12 o- The outer face is the marked face whose boundary has type 
(1,2,2). 



- 1625^^)5^3^ - 48(5^^)^5^^' - 2165^^) (.4^))^ - I20gi'\9i''r + 192(5^ 



and so on. 

The same treatment can be given to the other meromorphic parametrizations and obtain similar 
expressions : 

Mp) = -Vf(p+l)+T(^^+ gi'V 4.9« - 2,9^^ - 4.g,f + ... (5-21) 

.3(P) = ~Vt (2p + 1) + T ((3.g« + ) - Sgi'^ - 2gf^ - Agf) + ... (5-22) 

It is important to notice that these coefficients have to be found only once and thereafter they can 
be used to compute the T- expansions for the correlation functions without any further modification. 
Using these parametrizations and Eq. (5-15), we compute order by order the terms in the formal 
T-expansion for the disk amplitude with mixed boundary conditions i.e. 

W^i,2,3(.i,.2,.3) = E E (5-23) 

It is important to recall that 

00 ^ 

E (53)^" = ^lin^ ^ {Tr {M-^ M^^ Aq^)) . (5-24) 
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Here Trig's) denotes the contributions to the correlation function due to disks with v vertices and 
whose boundary conditions are given by an ordered sequence (ni,n2,?^3) of boundaries of types 1,2 
and 3 respectively. Since 172^'' — 52^'' the symmetry 

'{ni,7i2,n3)\93 '33 '53 J " ' (ng ,n2 ,"1 ) ^^3 '53 '53 ) [0-Zt,) 

holds. Moreover, the 7^(53)'s are given by power expansions in the couphngs of the cubic terms in 
the potentials, i.e 

T^i£) = E riM'^r'idPrHgi'Y^ (5-26) 

here has to be understood as the sum of the contributions to 7J containing rrii triangles of type 
i = 1,2, 3. It is worth seeing this explicitly for a couple of simple examples: 

7^1o,o)(5-^3) = - (4g^'' + 5^ + 25^'^) (5-27) 

7^1,2,0) (5l) = - {Wg^3^ + 45^ + 7^^ ') (5-28) 

The explicit enumeration of maps giving rise to the coefficients of these numbers is presented in fig|3] 
and fig|4j 

The contributions to W^i,2,3 up to three vertices (omitting terms that can be obtained by symmetry) 
are displayed in the following table. For clarity we use the simplified notation (jg^'' = gi, g^^^ = g2 and 

(3) 

53 = 53- 
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Boundaries \ Vertices 


1 


2 


3 


(1,0,0) 


1 


-4^3 - g2 - 2.g3 


-128.g? -64.9^52 -203153" 
-4.g| -96gfg3 -54515253 
-I65I53 -72515I -4O525I 
-645| 


(2,0,0) 





2 


645^ +245152 +45^ +4O5153 
+125253 +I65I 


(3,0,0) 








-325I -752 -1353 


(4,0,0) 








8 


(0,1,0) 


1 


-2.91 -92- 2^3 


-645? -405?5| -I6515I 
-45! -565^53 -42515253 
-I65I53 -56515I -4O525I 
-645| 


(0,2,0) 





1 


I65? +125152 +452 +I85153 
+125253 +165i 


(0,3,0) 








-751 -452 -753 


(0,4,0) 








2 


(1,1,0) 





1 


325? +175152 +452 +275153 
+125253 +165i 


(1,2,0) 








-IO51 -452 -753 


(1,3,0) 








2 


(2,1,0) 








-I651 -552 -953 


(3,1,0) 








4 


(2,2,0) 








1 


(1,0,1) 








-3252 -195152 -552 
-4I5153 -195253 -325I 


(2,0,1) 





-1 


I651 +752 +1453 


(3,0,1) 








-4 


(2,0,2) 








7 


(1,1,1) 








-251 -52 -453 


(2,1,1) 








1 


(1,2,1) 








2 



6 Conclusion and prospects 

In these notes, we found a recursive formula giving access to the so called complete disk amplitude 
in the multi-matrix model where an arbitrary number of matrices are coupled in chain. As in the 
two matrix model, this recursive formula only involves the spectral curve of the matrix model. Once 
the latter is known by the computation of the non-mixed disk amplitude, the recursive procedure 
presented in this paper explains how to get a disk with mixed boundary conditions by simple series 
expansion in terms of local variables on the spectral curve. This result proves a conjecture of Eynard 
[l3] and is the first step towards the generalization of the topological recursion formalism for the 
computation of mixed amplitude in the mutli-matrix setup. From the two matrix model experience 
|19j , it was also anticipated that the mixed correlation functions should follow from recursion relations 
taking into account all possible degenerations of the maps enumerated. 

This work obviously calls for generalizations and the computation of arbitrary mixed amplitude to 
any order in the expansion will probably be accessible by a recursive procedure similar to this 

one. This problem will be addressed in a forthcoming paper. 

On the other hand, a direct combinatorial understanding of these recursion relations is still missing 
and would deserve further investigations. Can we give any meaning to the so-called non-physical sheets 
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of the spectral curve? Can we get a universal factorization formula for arbitrary disk amplitude as in 
|17j ? These problems are definitely of high interest, not only for combinatorial motivations but also to 
investigate the generic structure underlying the mixed amplitudes in matrix models. This last point 
might be fundamental to understand the possible applications of these amplitudes to fields other than 
random matrix theories, such as topological string theories. 
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Appendix A Derivation of the loop equations 

A.l Loop equation for the complete mixed function 

Let us recall the derivation of the loop equations in the setup of I13j. 

1 1 ]_ 

X2 — M2 ■ ' ■ xj^ — Mj^ XX — M- 



Consider the change of variable 5Mi — — . . . — ztat ■ '^^'^ associated loop equation 



reads 



T 

(1-1) 

One can extract some interesting information from this equation. It can indeed read 



Res Wi^,,,^^{xi, . . . ,x^)V2{x2)dx2 = — Res 



Pi,i(xi, . . . ,x^f)V2{x2)dx2 



u Wi^3^4....,Af{xi,X3,...,Xj^)V2{x2)dx2 

+C1.2 Res — -— 

X2^oo Ci^2X2 — y (Xl) 

^ T Wi^...^.i{xi, . . . ,xu]Xi)V^{x2)dx2 

X2^oo Ci2X2—Y{xi) 

(1-2) 



One can remark that, for a function f{x) analytic in x — cxi 



Res ^/M^ ^Pol/(xi) (1-3) 

an-foo X — Xi xi 



which implies that 



Ci,2 Res Wi^,„^^r{xi,...,X^f)V2ix2)dX2 = Ci^2V2{x2)Wi^3A,...,Afixi,X3, . . . ,X^f) 

X2—>00 

T 

— PolWl.„.^f■l{xl, . . . , x^f■, xi)V2ix2) 

iV^ X2 



- Po\Pi,iixi,...,x^)V^{x2) 

X2 



Cl,2X2=Y(xi) 



(1-4) 

Using the notations of the preceding section, one gets 

Ci,2 Res Wi^,„x{xi,...,X^f)V2{x2)dX2 = -Pi^2{xi,X2,X3,...,X^) 

-ci,2W3,4,...x(a;3, ■ ■ • ,a;7v) 

+Ci,2V2'(X2)VFi,3.4,....AA(a;i,a;3, . . .,xj^) 

T 

+ J^P2,2-l{.Xi,X2,X3, . . .,X^f). 



(1-5) 



A. 2 Hierarchy of loop equations 



Let us now consider the change of variable 5M2 = x^-m^ x^-au ' ' ' xi-Mi ■ ^^'^ corresponding loop 
equation is 

- l^/Tr ^—V^{M2) ^\+ci,2W^3,4,...,AAfe,-.-,a;A^) 

iV \ xi- Ml X3 - A/3 x^f - A'/v / 

+C2,3W^i,4,...,AA - (ci,22:i + 02,32^3) W^i,3,4,.. . .A/" (a;i , a;3 , . . . , x^r). 
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(1-6) 



The first term is actually 



-T^--- -Tr~/^ Res Wi^...^m{xi,- ■ ■ ,xM)V2{x2)dx2 (1-7) 



which implies, using Eq. ( 1-5 ) 



{X3 - Xsixi)) Wi^3^...^^f{xi,X3, . . ■,X^f) 



where 



Wi^4^,„Xixi,X4, ■ ■ ■,XJ^) 

Pol Wy,...,M\Xx,-..,X}4) 

X1.X2 C1.2C1.3 

- Pol W^1,...,AA(xi, . . . ,X7V) 



XI, X2 C2.3 
iV^ X2 Ci,2C2,3 



VFi,...x(xi, . . ■ ,X7V) 



X3 



(1-8) 

^2(^2) - Ci,2Xi 



C2,3 



(1-9) 



Let us now proceed by induction. Let < fc < A/" and assume that one has 



{Xj - Xj{xi))Wi^j^,„^j^{xi,Xj, . . . ,Xj^) = 

= M^ij+i,...,A"(xi,Xj+i, . . . ,XAr) 



Pol /lj-l(xi, . . . , Xj_i)H/"i,...,A/-(xi, . . . ,Xa^) 



Xi,...,Xj-X 

T 



{xi=XiYi^l 



Pol /2j-i(x2, . . . ,Xj-l)Wl„„^^f{xl, . . .,x^) 



X2,...,Xj-l 



(1 - 10) 



{xi—Xi}-'^^^ 



for j < k. We now show that it is also true for j = fc + 1. 
This property for j = k implies that 



(1-11) 



Res Vl{xk)dxkWi^k,...,Mi^i^^k, x^) 



Res 

Xk 



Res 

Xk 



Vl{xk)dxkWi^k+i,...Mixi,Xk+i, ■ ■ ■ ,x^f) 

{Xk - Xfc(xi)) 

Vl{Xk)dXk [PoUi,...,2;fc_i /l,fc-l(xi, . . ■,Xk-l)Wi^,,,^^r{xi,.. . ,Xf^)]^^_^^^k-l 



{xk - Xk{xi)) 

T Vl{xk)dxk [Polx2....,Xk-i /2,fc-i(x2, . . . ,Xk~i)Wi^...x{xi, . . .,x^f)]. _ .k^i 

Res ^ '~ '-^'=^ 

N'^ Xk {xk-Xk{xi)) 



Res Vk{xk)dxkWi^k.....Af{xi,Xk, ■ ■ . ,x^f) = 

Xk 

Vk{xk)Wi^k+i,...Mixi,Xk+i, ■ ■ ■,x^f) 

Pol Vj;(xfc)/i,fc_i(a;i, . . . ,Xk-l)Wl^,„^^f{xl, . . . ,xaa) 

Xi,...,Xk 



T 



Pol Vfe(Xfc)/2.fe-l(x2, ■ • • ,Xfc_i)Wi....,AA(xi, 

X2,....Xk 



{Xi—X 

Xat) 



1* 

li = 2 
{xi=Xi}'^^2 
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(1-12) 



Let us consider the induction hypothesis for j = fc — 1, it gives 
Res dxk-iXk-iWi k-i,...,Ar{xi,Xk-i, ■ • ■,xj^) = 

Xk-l 

Xk-iWi^k.....xixi,Xk, . . . , x^r) 

Pol Xk-lfl,k-2{xi, . ■ ■ ,Xk-l)Wi^,„J^{xi, . . .,XJ^) 

Xi,...,Xk-l 



T 



Pol Xk-lf2,k-2{xi, . . . , Xk-l)Wi^...^{xi, Xaa) 

X2,...,Xk-l 



{xi=Xi}- 



(1-13) 

The loop equation corresponding to SMk = ^,^^,^^+1 • • • xm-m^t ^T^mT ''^^^^ (^^'^ k < Af) 



= ^ / Tr ^-^VliMk) 



1 



1 



N \ xi-Mi 'xk+i - Mk+i XM-M^ ^ 

-Ck-i,k R-es dxkdxk-iXk-\Wik-i,...j^{xi,Xk-i,...-,xj^) 

Xk-l,Xk 

+Ck,k+lWi^k+2,...,Ar{xi,Xk+2, ■ ■ ■,XX) 
-Ck,k+lXk+lWi^k+l,...,M{xl,Xk+l, . . ■,Xj^). 

(1-14) 

Plugging in the expression found earlier, one gets 

{Xk+l - Xk+l{Xl)) Wi^k+l,...,M{Xi,Xk+l, ■ ■ ■,x^r) = 

Wi^k+2,...,Ar{xi,Xk+i,...,x^f) 

Pol f\,k{xi,. ■ ■ ,Xk-l)Wi^,„^J^{xi, . . .,XJ^) 

Xl,...,Xk-l 



+ 



{xi=Xi}'?^ 



Pol f2,k{X2, Xk-l)Wi^...x{xi, Xm) 

X2,--.,Xk-l 



(1-15) 



{Xi—Xi}'^_ 



with 



and 



Xk+l 



Ki^k) - Ck-l,kXk-l 



Ck,k+1 



, / N K'(a;ft)/i fe_i — Ck-i,kXk-iXkfi,k-2 

fl,k{Xl, Xk-l) = — — ■ ■ 

Cfc,fe+1 

J. . ^ ykiXk)f2,k-l - Ck-l.kXk-lXkf2,k-2 
f2,k{xi,...,Xk-{) = — — ■ . 

Cfc,fc+1 



Together with the initial equation for Wi^{xi,x_\f), this ends the proof. 



(1-16) 
(1-17) 

(1-18) 
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